Monte Carlo tomographic reconstruction in SPECT impact of bootstrapping and number of generated events

In Single Photon Emission Computed Tomography (SPECT), 3D images usually reconstructed by performing a set of bidimensional (2D) analytical or iterative reconstructions can also be reconstructed using an iterative reconstruction algorithm involving a 3D projector. Accurate Monte Carlo (MC) simulations modeling all the physical effects that affect the imaging process can be used to estimate this projector. However, the accuracy of the projector is affected by the stochastic nature of MC simulations. In this paper, we study the accuracy of the reconstructed images with respect to the number of simulated histories used to estimate the MC projector. Furthermore, we study the impact of applying the bootstrapping technique when estimating the projectorComment: 15 pages, 9 figures, 2 table

Fermion-scalar interactions with domain wall fermions

Domain wall fermions are defined on a lattice with an extra direction the size of which controls the chiral properties of the theory. When gauge fields are coupled to domain wall fermions the extra direction is treated as an internal flavor space. Here it is found that this is not the case for scalar fields. Instead, the interaction takes place only along the link that connects the boundaries of the extra direction. This reveals a richness in the way different spin particles are coupled to domain wall fermions. As an application, 4-Fermi models are studied using large N techniques and the results are supported by numerical simulations with N=2. It is found that the chiral properties of domain wall fermions in these models are good across a large range of couplings and that a phase with parity-flavor broken symmetry can develop for negative bare masses if the number of sites along the extra direction is finite.Comment: LaTeX, 17 pages, 8 eps figures; comment regarding the width of Aoki phase added in sec. 3; references adde

Effect of noise and modeling errors on the reliability of fully 3D Monte Carlo reconstruction in SPECT

We recently demonstrated the value of reconstructing SPECT data with fully 3D Monte Carlo reconstruction (F3DMC), in terms of spatial resolution and quantification. This was shown on a small cubic phantom (64 projections 10 x 10) in some idealistic configurations. The goals of the present study were to assess the effect of noise and modeling errors on the reliability of F3DMC, to propose and evaluate strategies for reducing the noise in the projector, and to demonstrate the feasibility of F3DMC for a dataset with realistic dimensions. A small cubic phantom and a realistic Jaszczak phantom dataset were considered. Projections and projectors for both phantoms were calculated using the Monte Carlo simulation code GATE. Projectors with different statistics were considered and two methods for reducing noise in the projector were investigated: one based on principal component analysis (PCA) and the other consisting in setting small probability values to zero. Energy and spatial shifts in projection sampling with respect to projector sampling were also introduced to test F3DMC in realistic conditions. Experiments with the cubic phantom showed the importance of using simulations with high statistics for calculating the projector, and the value of filtering the projector using a PCA approach. F3DMC was shown to be robust with respect to energy shift and small spatial sampling off-set between the projector and the projections. Images of the Jaszczak phantom were successfully reconstructed and also showed promising results in terms of spatial resolution recovery and quantitative accuracy in small structures. It is concluded that the promising results of F3DMC hold on realistic data set

Reconstruction tri-dimensionnelle complete d'images en spect-ct par modelisation Monte-Carlo

présenté par Z. El Bitar, proceedings sous forme de CDEn tomographie d'émission monophotonique (SPECT), les images 3D normalement reconstruites par des algorithmes de reconstruction analytiques ou itératives bidimensionnelles (2D) pourraient aussi bien être reconstruites avec des algorithmes de reconstruction itérative (3D) qui permettent de compenser les effets physiques perturbant le processus de formation de l'image notamment l'atténuation et la diffusion. Nous avons étudié une technique de reconstruction 3D complète (F3DMC) (Lazaro et al. NIM 2004), dans laquelle le projecteur 3D impliqué dans la reconstruction est estimé par des simulations Monte-Carlo effectuées à partir de données tomodensitométriques du patient

Targeted fully 3D Monte Carlo reconstruction in SPECT

PCSV, présenté par Z. El Bitar, soumis aux proceedingsFully 3D Monte-Carlo (F3DMC) reconstruction consists in calculating a fully 3D object-specific system matrix using Monte-Carlo simulations and inverting it using an iterative approach. To reduce the large amount of disk space required by this approach, we derived a targeted F3DMC approach (TF3DMC) in which the volume to be reconstructed is irregularly sampled, so that pre-identified functional regions of interest are reconstructed using fine sampling while regions with non-specific activity or without any particular interest are coarsely sampled. This method was assessed using simulated and real SPECT data of a phantom filled with Tc99m. The GATE Monte-Carlo simulator was considered to simulate the phantom data and to calculate the system matrices needed for the reconstruction of the simulated and of the real SPECT data. Activity ratios measured in TF3DMC images were compared with those measured on F3DMC and OSEM images corrected for scatter, attenuation and detector response function. TF3DLMC yielded errors less than 10% in activity ratio estimates in hot regions, while errors with quantitative OSEM were between -21% and -3%. The space needed to store the system matrix was divided by a factor from 3.5 to 9.4 compared to F3DMC, for similar or even better accuracy in activity ratio estimates. These results suggest that TF3DMC can be made practical and outperforms F3DMC and OSEM in terms of quantitative accuracy

Monte Carlo Simulation With The GATE Software Using Grid Computing

DémonstrationInternational audienceMonte Carlo simulations needing many replicates to obtain good statistical results can be easily executed in parallel using the "Multiple Replications In Parallel" approach. However, several precautions have to be taken in the generation of the parallel streams of pseudo-random numbers. In this paper, we present the distribution of Monte Carlo simulations performed with the GATE software using local clusters and grid computing. We obtained very convincing results with this large medical application, thanks to the EGEE Grid (Enabling Grid for E-sciencE), achieving in one week computations that could have taken more than 3 years of processing on a single computer. This work has been achieved thanks to a generic object-oriented toolbox called DistMe which we designed to automate this kind of parallelization for Monte Carlo simulations. This toolbox, written in Java is freely available on SourceForge and helped to ensure a rigorous distribution of pseudo-random number streams. It is based on the use of a documented XML format for random numbers generators statuses

How Can Autonomous Road Vehicles Coexist with Human-Driven Vehicles? An Evolutionary-Game-Theoretic Perspective

The advent of highly automated vehicles in the form of autonomous road vehicles (ARVs) is bound to bring about a paradigm shift in road user interaction, especially that between ARVs and human-driven vehicles (HDVs). Previous literature on the game-theoretic interaction between ARVs and HDVs tends to focus on working out the best possible strategy for a single interaction, i.e. the Nash equilibrium. This position paper sets out to demonstrate the importance and potential impact of applying evolutionary game theoretic principles to what is effectively a dynamic population driven by evolutionary forces – the population of road users. We demonstrate using theoretical scenarios that simply maintaining Nash equilibria does not guarantee evolutionary success. Instead, ARVs must enjoy a demonstrable advantage over other road users when few in numbers. Otherwise, their uptake will slow down and eventually reverse. We argue that the same selection factors which influence the success of living populations in the natural world also influence the success of the different vehicle types and driving styles in the road user population, including ARVs. We demonstrate this by assigning an expected fitness score to each vehicle in a one-to-one interaction, such as at a junction. This fitness score is dependent on driver, rider and economic costs incurred by the vehicle and/or its occupant(s) during interaction. In turn we show that ARV and transport system designers need to ensure that the fitness score of their systems create evolutionary stability

Two-dimensional Lattice Gross-Neveu Model with Wilson Fermion Action at Finite Temperature and Chemical Potential

We investigate the phase structure of the two-dimensional lattice Gross-Neveu model formulated with the Wilson fermion action to leading order of 1/N expansion. Structural change of the parity-broken phase under the influence of finite temperature and chemical potential is studied. The connection between the lattice phase structure and the chiral phase transition of the continuum theory is clarified.Comment: 42 pages, 20 EPS figures, using REVTe

Dynamical Wilson fermions and the problem of the chiral limit in compact lattice QED

We compare the approach to the chiral transition line ~\kappa_c(\bt)~ in quenched and full compact lattice QED with Wilson fermions within the confinement phase, especially in the pseudoscalar sector of the theory. We show that in the strong coupling limit ($\beta =0$) the quenched theory is a good approximation to the full one, in contrast to the case of $\beta =0.8$. At the larger $\beta$-value the transition in the full theory is inconsistent with the zero--mass limit of the pseudoscalar particle, thus prohibiting the definition of a chiral limit.Comment: 13 pages LaTeX (epsf), all figures include

Complete Nucleotide Sequence of Plasmids of Two Escherichia coli Strains Carrying blaNDM–5 and blaNDM–5 and blaOXA–181 From the Same Patient

Aim of this study was to genetically characterize two carbapenemase-producing Escherichia coli strains obtained from a pediatric patient affected by diarrhea, expressing OXA-181 and/or NDM-5 type enzymes. The above microorganisms were collected in the same Desenzano hospital (Northern Italy) where the blaNDM–5 gene was detected for the first time in Italy 3 years ago. One strain (5P), belonged to sequence type ST405/ST477 (according to Pasture/Oxford schemes) and serotype O102:H6. It was characterized by a 130562 bp multi-replicon plasmid IncFII/IncFIA/IncFIB (pVSI_NDM-5) enclosing two main antibiotic resistance islands: (i) ARI-I, 10030 bp in size, carried genes coding for β-lactam- (blaOXA–1, blaCTX–M–15), fluoroquinolone/aminoglycoside- (aac(6′)-lb-cr) and phenicol- resistance (catB3), (ii) ARI-II, 15326 bp in size, carried genes coding for sulfonamide- (sul1), β-lactam- (blaNDM–5, blaTEM–1B), phenicol- (catB3), trimethoprim- (dfrA17), antiseptic- (qacEΔ1), and aminoglycoside- (aadA5, rmtB) resistance. The other isolate (5M), belonged to sequence type ST2659/ST759 and serotype O50/02:H18, and carried four plasmids: a 153866 bp multi-replicon IncFII/IncFIA/IncFIB (pISV_IncFII_NDM-5), an 89866 bp IncI1 plasmid, a 51480 bp IncX3 plasmid (pISV_IncX3_OXA181), and a 41143 bp IncI plasmid (pISV_IncI_CMY-42). pISV_IncFII_NDM-5 carried two main antibiotic resistance islands: (i) ARI-III, 12220 bp in size, carried genes coding for β-lactam- (blaOXA–1), fluoroquinolone/aminoglycoside- (aac(6′)-lb-cr), tetracycline- (tet(B)) and phenicol- resistance (catB3, catA1), and ii) ARI-IV, 26527 bp in size, carried determinants coding for macrolide- (erm(B), mph(A)), sulfonamide- (sul1), beta-lactam- (blaNDM–5, blaTEM–1B), trimethoprim- (dfrA14, dfrA12), antiseptic- (qacEΔ1), and aminoglycoside- resistance (aadA5). pISV_IncI_CMY-42 harbored the blaCMY–42 gene coding for beta-lactam resistance, pISV_IncX3_OXA181 harbored genes encoding fluoroquinolone- (qnrS1) and beta-lactams- resistance (blaOXA–181). In conclusion, the detection of two different NDM-5 E. coli strains from a pediatric patient with a history of travel to the Far East countries strongly highlight an increasing trend and risk of importation from such areas